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Abstract 

I The chiral Lagrangian for the Symanzik action through 0{a'^) for baryons is obtained. We 

' consider two flavor unquenched and partially quenched lattice theories, allowing for mixed actions 

in the latter. As an application, we calculate masses to O(a^) for the nucleons and deltas, and 
pi |. investigate the corrections due to the violation of 0(4) rotational invariance. These results are 

, contrasted with those in the meson sector for lattice simulations using mixed and unmixed actions 

of Wilson and Ginsparg- Wilson quarks. 
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I. INTRODUCTION 



Lattice gauge theory provides first principles calculations of strong interaction physics, 
where QCD is non-perturbative, and quarks and gluons are confined in color- neutral 
hadronic states. These calculations, however, are severely limited by available computing 
power, necessitating the use of quark masses that are much larger than those in reality. 
To make physical predictions, one must extrapolate from the quark masses used on the lat- 
tice to those of nature. A model independent tool for this extrapolation is to study QCD 
at hadronic scales using its low-energy effective theory, chiral perturbation theory (^PT). 
Because xPT provides a systematic expansion involving niq/ Kqcd-, one can understand how 
QCD observables behave, in principle, as functions of the quark mass. To address the 
quenched and partially quenched approximations employed by lattice calculations, xPT has 
been extended to quenched chiral perturbation theory (OyPT) Ei IE 1^ ^-iid partially 
quenched chiral perturbation theory (PQxPT) jsi 0, H, 0, El • 

While lattice calculations are limited to unphysically large quark masses, they are re- 
stricted further by two additional parameters: the size of the lattice L, that is not con- 
siderably larger than the system under investigation; and the lattice spacing a, that is not 
considerably smaller than the relevant hadronic distance scale. To address the issue of finite 
lattice spacing, xPT has been extended (following the earlier work of 0, [ ^ in the meson 
sector to 0{a) for the Wilson action Il4l | , a nd for mixed lattice actions 15[. Corrections 
at 0{a?') in xPT have been pursued |l6l |l7|. There have also been parallel developments 
in addressing lattice spacing artifacts in staggered yPT for mesons [3, [3, HO], and heavy 
mesons j2l|; and in twisted mass QCD [j^ liil E3 . lisl . li^ . Corrections to baryon observ- 
ables in xPT and PQxPT, too, have been recently investigated 27, 2^. To consider finite 
lattice spacing corrections, one must formulate the underlying lattice theory and match the 
new operators that appear onto those in the chiral effective theory. This can be don e b'v 
utilizing a dual expansion in quark mass and lattice spacing. For an overview, see 
Following 3, 131, we assume a hierarchy of energy scales 



mq < Kqcd < ^. (1) 

We shall further choose a power counting scheme in which the small dimensionless expansion 
parameters are^ 

Thus we have a systematic way to calculate a-dependent corrections in xPT for the observ- 
ables of interest. Such expressions allow one to perform the quark mass extrapolation before 
the continuum extrapolation. 

In this work we address the extension of xPT and PQxPT at finite lattice spacing to 
O(a^) in the baryon sector. Specifically we detail the operators which must be included to 
determine the baryon masses to 0{e'^) in the above power counting. Unlike the extension of 
XPT and PQxPT to 0{a^) in the meson sector, the baryon sector suffers from a proliferation 
of new operators. Despite this fact, the number of free independent parameters entering 



^ In practice the power counting scheme and subsequent ordering of the dual expansion in quark mass and 
lattice spacing should be organized based the on the actual sizes of niq and a. 
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expressions for the masses is still relatively small. Moreover in contrast to the meson sector 
at 0(0^), heavy baryon operators that break the 0(4) rotational symmetry of Euclidean 
space are required at this order in the chiral expansion. Such operators are required for 
particles that are heavy compared to Aqqd- 

This paper has the following organization. First in Sec. Ill A[ we review the Symanzik 
Lagrangian at O(a^) for a partially quenched, mixed lattice action, where the valence and 
sea quarks are either Wilson or Ginsparg- Wilson fermions. We focus on the symmetries of 
the Symanzik Lagrangian because these are essential in constructing the effective theory. 
Next in Sec. Ill Bl we review the finite lattice spacing corrections at 0{a) in the meson sector 
of PQxPT. Higher-order corrections are not needed for baryon observables to the order we 
work. In Sec. Ill Cl we extend heavy baryon PQxPT to O(e^). This includes the addition of 
operators at O^mga) and 0{a^). Apphcations of this development are pursued in Sec. IIIII, 
where we obtain the lattice spacing corrections to the masses of nucleons and deltas. The 
unquenched two-flavor theory is addressed in Appendix^ Corrections from 0(4) breaking 
operators are treated in detail for particles of spin less than two in Appendix^ A summary 
(Sec. I1V|) highlights the lattice spacing corrections in the baryon sector. Here we contrast 
the number of independent parameters entering expressions at O(e^) for baryon masses in 
the various mixed and unmixed partially quenched theories, as well as in the unquenched 
theory. 

II. PQxPT AT 0(0^) 

To extend the baryon chiral Lagrangian to O(a^), we flrst review the Symanzik La- 
grangian at O(a^) and the construction of PQxPT in the meson sector. In this and the 
following sections, we consider a partially quenched theory with a mixed action. The result 
for unquenched simulations is contained in Appendix 1X1 

A. Symanzik Lagrangian 

The Symanzik action is the continuum effective theory of the lattice action [s^ |^ . As 
such it is constructed from continuum operators based on the symmetries of the underlying 
lattice theory. The Symanzik Lagrangian is organized in powers of the lattice spacing a, 
namely 

C = C^^^+aC^^^ + a^C^^^ + ..., (3) 

where C^""^ represents contributions for dimension-ra operators.^ In the continuum limit, 
a — i> 0, only the dimension-four operators survive. The Symanzik Lagrangian for the Wilson 
action was discussed to C(a^) in ji^l and the analysis to 0{a) was refined by jl^. Here 
we consider the general case of a mixed action in partially quenched QCD (PQQCD). Such 
an action allows the valence and sea quarks to have different masses. Additionally the va- 
lence and sea quarks can be different types of lattice fermions. The mixed lattice action 
has the usual parity invariance, charge conjugation invariance and SU (iVc) gauge symmetry. 



^ Not all of the a-dependence is parametrized in Eq. (|3Jl. The coefficients of terms in the Lagrangian /I'^") 
depend upon the gauge coupling and thus can have a weak logarithmic dependence on a. Such dependence 
is beyond the scope of this work. 
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Accordingly the Symanzik Lagrangian L respects these symmetries order-by-order in a. Be- 
cause spacetime has been discretized, the 0(4) rotational symmetry of continuum Euclidean 
field theory has been reduced to the hypercubic group. 

The fiavor symmetry group of the mixed lattice action G is generally a direct product 
of fiavor symmetry groups for the particular species of fermion at hand The mixed 

action contains different forms of the Dirac operator for each species of fermion, thus there 
is no symmetry transformation between the valence and sea sectors. At zero quark mass, 
the fiavor symmetry group of the mixed lattice action is 

G = G valence ® Gsea., (4) 

where for two quark fiavors 

Gvalence = SU{2\2)l ® SU{2\2)r, (5) 

and 

aea = SU{2)l ® SU{2)r. (6) 

For simplicity we have included the ghost quarks in the valence sector. The quark mass 
term breaks these respective chiral symmetries in the valence and sea sectors down to vector 
symmetries. 

For Wilson fermions [s^l, the breaking of chiral symmetry persists even at zero quark 
mass due to lattice discretization effects. These effects can be tamed to enter at O(a^) by 
the so-called 0{a) improvement. For fermions satisfying the Ginsparg- Wilson relation [35| 
(e.g., Kaplan fermions or overlap fermions js^), the chiral symmetry on the lattice at 
zero quark mass remains exact at finite a jHI. For a recent review of chiral symmetry in 
lattice theories, see !39]. We have left unspecified which species lives in which sector of the 
theory. As a result of this generality, we are building the Symanzik Lagrangian for the four 
possible combinations of valence and sea fermions. Part of this generality is an academic 
pursuit because we do not anticipate lattice calculations employing Wilson valence quarks 
in a Ginsparg- Wilson sea, whereas GW valence quarks with dynamical Wilson fermions is 
a scenario with potential computational benefits. It is likely, however, that Wilson quark 
masses may never reach the chiral regime ji^] in which case one must deal with alternate 
methods to modify 41 1 or to improve '421] baryon xPT. These reservations aside, our goal is 



to contrast the situation at O(a^) in the baryon sector with that of the mesons. We will find 
comparatively that the mesons are rather special with respect to lattice spacing corrections 
in xPT. 

Having discussed the symmetries of C, we now turn to the contributions near the contin- 
uum limit. The terms of the C^^^ Lagrangian are the most familiar and so we discuss them 
first. At dimension four, we have the familiar kinetic and Dirac mass terms for the quarks. 
To be specific, we work in a partially quenched theory for two light fiavors. The quark part 
of the leading-order Symanzik Lagrangian reads 

C^"^^ =Q^Q + QmQQ. (7) 

The six quarks of PQQCD are in the fundamental representation of the graded group 
and appear in the vector 

Q = {u,d,i,l,u,df , (8) 
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which obeys the graded equal-time commutation relation 

Qr(x)gf (y) - i-iy'"'^Qfiy)Qn^) = s^n^.s^-^ - y), (9) 

where a and (5 are spin, and i and j are flavor indices. The vanishing graded equal-time 
commutation relations can be written analogously. The grading factor 

/ 1 for A; = 1,2,3,4 

I for A; = 5,6 ' (1°) 

incorporates the different statistics of the quarks. The quark mass matrix in the isospin 
limit {md = mu) of S't/(4|2) is given by 

niQ = dia:g{mu,mu,mj,mj,mu,mu). (11) 

In the limit a — 0, and when rrij = rriu, one recovers the isospin limit of QCD. The symmetry 
of C^^^ in the zero mass limit is SU{4:\2)l ® SU{4:\2)r, i.e. the effects of the mixed action do 
not show up in the effective theory at leading order. The mass term in C^^^ explicitly breaks 
the graded chiral symmetry down to the graded vector symmetry SU{4\2)y. 

Let us next consider the dimension-flve Lagrangian. After fleld redeflnitions, it consists 
of only the Pauli term given by 

jC^'^ = cswQ(Tf.uG^,WQQ. (12) 

This term breaks chiral symmetry in precisely the same way as the quark mass term. The 
Sheikholeslami-Wohlert (SW) !32| coefficient is csw and it is accompanied by the Wilson 
matrix wq defined by 

wq = diag(w^„ Wy, Ws, Ws, Wy, Wy). (13) 

This matrix accounts for the chiral symmetry properties of the mixed action. If the quark 
Qi is a Wilson fermion, then (wq)^ = 1. Alternately, if Qi is of the Ginsparg- Wilson (GW) 
variety then (wq)^ = 0. Since one expects simulations to be performed with valence quarks 
that are all of the same species as well as sea quarks all of the same species, we have labeled 
the entries in Eq. (jl3p by valence (v) and sea (s) instead of fiavor.^ 

To work at (9(a^), we must consider the dimension-six Lagrangian. The exact form of C^^^ 
is irrelevant in constructing the chiral effective theory. We need only know which symmetries 
are broken and how. We explicitly list terms of the Lagrangian when it is illustrative to 
do so. In describing these operators, we introduce the flavor matrix wq = 1 — wq, which 
projects onto the GW sector of the theory. 

The terms of C^^^ fall into flve classes. The flrst class of operators consists of higher- 
dimensional quark bilinears. These operators do not break chiral symmetry but have the 
flavor symmetry of the mixed action, i.e. there are distinct operators involving only the 
valence sector and only the sea sector. Typical examples are Q Ip^wq Q, and Q Ip^wqQ. 



^ It is conceivable that three flavor simulations might employ Wilson up and down quarks but with a GW 
strange quark in order to eliminate the rather large 0{ams) corrections. The development in this work 
for partially quenched theories is general enough to handle this situation by a modification of Eq. I 
thank W. Detmold for bringing this to my attention. 
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TABLE I: Description of the classes of operators in C^^\ Each representative class of operator is 
classified as to whether it is bilinear, chiral symmetry breaking (xSB), or 0(4) breaking. Operators 
m £(6) that are not bilinears are thus four-quark operators. Additionally the flavor matrices 
involved for the mixed action are listed. 



Class 


Bilinear? 


XSB? 


0(4) breaking? 


Flavor Structure 


1 


Yes 


No 


No 


WQ, WQ 


2 


Yes 


Yes 


No 


mQ®WQ, ruQ <^wq 


3 


No 


No 


No Wq 


^ WQ, WQ (g) WQ, WQ (gl WQ 


4 


No 


Yes 


No 


WQ (g) WQ 


5 


Yes 


No 


Yes 


WQ, WQ 



In class two, there are quark bilinear operators that break chiral symmetry. Simple di- 
mensional counting indicates that there must be an ruQ insertion for these operators to be 
in C^^\ Examples include Q niQ D'^wq Q, and Q tuq D'^wq Q. Class three consists of all 
four-quark operators that do not break the chiral symmetry of the valence and sea sec- 
tors. These class three operators come in three forms due to the different sectors of the 
theory, e.g. {QwQ-i^QY, {Qwq^i^QY, and {QwQ'yf,Q){QwQ'y^Q). Class four operators 
are four-quark operators that break chiral symmetry. Since there are four fields, chiral 
symmetry is broken in a different way than the quark mass term. A typical class four 
operator is {Q wqQY . Unlike class three, there are no class four operators involving the 
flavor matrix wq. Finally class five operators are those that break 0(4) rotation symmetry. 
There are two such operators for the mixed lattice action, namely Q wq j^D^D^D^Q, and 
QujQ'-y^D^D^D^Q. These classes are summarized in Table HI 



B. Mesons 



For massless quarks at zero lattice spacing, the Lagrangian in Eq. (|7j) exhibits a graded 
symmetry SU{4:\2)l S'[/(4|2)r ® U{l)v that is assumed to be spontaneously broken to 
SU (4|2)y (g) U {l)v- The low-energy effective theory of PQQCD that results from perturbing 
about the physical vacuum is PQxPT. The pseudo-Goldstone mesons can be described at 
0{e^) by a La grangian that accounts for the two sources of explicit chiral symmetry breaking: 



ral syrnr 
313: 



the quark mass term in Eq. (|7j), and the Pauli term in Eq. ()12|) ll2L ll 

C = ^str [df.E'^ df^Tj) - str ^mgS''" + m^S j - aXa str (^wqTj'' + WqS j (14) 

where 

E = exp(H^)=^^ (15) 

/ = 132 MeV, and the str() denotes a graded flavor trace. The M, M, and \ are matrices of 



Y are 

pseudo-Goldstone bosons and pseudo-Goldstone fermions, see, for example, |46|. Expanding 
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the Lagrangian in ()14p one finds that to lowest order mesons with quark content QQ' have 

4 



mass^ 



QQ' = -p \^rn{mQ + niQi) + a\a{wQ + Wq.)] . (17) 

The flavor singlet fleld is rendered heavy by the U{1)a anomaly and has been integrated out 
in PQ;yPT, however, the propagator of the flavor-neutral fleld deviates from a simple pole 
form lnJl. For a,b = u, d,j, /, u, d, the leading-order rjarji, propagator is given by 

^ iea6ab _ i_ egtb [g^ - ni]^) 

g2 _ + 2 (g2 _ m2^^ + ie) {f - ml + ze) ' ^ ^ 

where 

e, = (-l)l+^^ (19) 
The flavor neutral propagator can be conveniently rewritten as 



^VaVb ~ ^aSabPa + ^a^h^ab (Pa, Ph) , (20) 



where 



g2 — ml + ie' — ml + is ' 



Uab {A,B) 



Va J] A Vb Jj R 

— mf.^ mf.^ — mf. 

Va 'lb Vb 'la 



(21) 



At 0{e^), one has contributions to the meson Lagrangian from operators of 0(p^), 
0{p'^ mq), and 0{mQ). These are the Gasser-Leutwyler terms. Additionally there are terms 
of 0{p'^a), O^mqa), and C(a^) that are generalizations of the Gasser-Leutwyler terms for 
the Symanzik lattice action. These have been determined in and we do not duplicate 

them here. Since our concern does not lie in the meson sector, mesons will only enter via 
loop calculations. Retaining the meson masses to O(e^) in a typical baryon calculation leads 
to corrections of 0{e^) or higher. These are beyond the order we work, thus the Lagrangian 
in Eq. (fT^ is sufficient for our purposes. 



C. Baryons 

Having reviewed the Symanzik Lagrangian through O(a^) and the relevant pieces of 
meson PQxPT at finite a, we now extend PQxPT in the baryon sector to O(a^). First let 
us detail the situation at 0{a). In S'f/(4|2) PQxPT, the spin-| baryons are embedded in the 



The quark masses niQ above are not those customarily used on the lattice |12L llTj , because one usually 
defines the quark mass in terms of a critical parameter for which the meson masses vanish. If this is 
indeed the way one defines the renormalized quark mass, then the parameter Aa in the Lagrangian can 
be set to zero and the issue of a-dependent loop-meson masses in baryonic observables will never confront 
us. Allowing for other definitions of the lattice renormalized quark mass to avoid this fine-tuning, we keep 
Aa ^ and deal with the possibility that the loop- meson masses remain a-dependent. Notice this issue 
only arises for Wilson quarks. 
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70- dimensional super-multiplet B^^^, that contains the nucleons, while the spin-| baryons 
are embedded in the 44-dimensional super-multiplet T^^'', that contains the deltas 0, 0|- 
To the free Lagrangian for the B^^'' and T*-^'^ fields is given by |3. li?! li?! 



C = i{Bv VB) - 2a M (BBM+) - 2pM (BM+B) - 2a m (BB) str {M+) 
-2aw {BBW+) - 2pw (BW+B) - 2aw (BB) str 

+% {T^v ■ VT^) + A (T^%) + 27m (T^M+T^) - 2aM (T^T^) str {M+) 
+2^w - 2aw (T^T^) str (W+) , (22) 



where the mass operator is defined by 



M± = l {^^rriQe ± ^niQ^) , (23) 



and the Wilson operator is defined by^ 

Wi = i^WQ^^ ± ^wqC) . (24) 



Here A ~ e is the mass splitting between the 70 and 44 in the chiral limit. The parenthesis 
notation used in Eq. (j^^ is that of |J| and is defined so that the contractions of flavor 
indices maintain proper transformations under chiral rotations. Notice that the presence 
of the chiral symmetry breaking SW operator in Eq. ((7j) has lead to new 0{a) operators 
and new dimensionless constants aw, Pw, <^w, Iw, and aw in Eq- (f^ - The Lagrangian 
describing the interactions of the B^^'^ and T^^'' with the pseudo-Goldstone mesons is 



£ = 2a {BS^BA^)+2(3 {BS^A^B)-2n {T,S^A^%) + [{T.A^B) + {BA,%)] . 



(25) 



The axial- vector and vector meson fields A^ and are defined by: A^ = ^ {^d^C,^ — ^^9^^) 
and = \ (^9^^^ + ^^<9^^)- The latter appears in Eq. for the covariant derivatives of 
Bijk and Tijk that both have the form 



{V,B),,u = d,B,,k + {y^hiBi.k + {~y'^''''-''-\V,),mBimk + (-)^'''+''^-^^'"^+''"^(V^)fcni3,,„. (26) 



The vector S*^ is the covariant spin operator |48|, |49|, |50[. The interaction Lagrangian in 
Eq. also receives finite a corrections. In calculating the octet and decuplet masses, 
however, these lead to effects that are of 0{e^) or higher. 

At 0{e'^), there are contributions to the PQxPT Lagrangian from two insertions of the 
mass operator + , and contributions from two insertions of the axial current A^. The 
former contribute to the baryon masses at tree level, the latter at one-loop level. These 
operators have been written down in 51, 5^ 5^. There are also operators with an insertion 



^ Technically the SW coefficient csw: with its possible weak logarithmic dependence on a, also enters into 
the spurion construction. We are ignoring this dependence and do not distinguish between spurious of 
the same form which come from different operators in the Symanzik action largely because at O(a^) the 
analysis becomes extremely cumbersome. 
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of V ■ A® A4^. These do not contribute to the masses at C>(£^). To extend baryon PQxPT 
for mixed lattice actions to O(e^), we must first include all higher-order operators that 
are linear in a. Secondly we must map the operators in C^^^ onto the baryon sector. To 
achieve the former, we realize that the relevant operators can be formed from the existing 
0{a) operators by insertion of the mass operator (insertion of a derivative is ruled 
out because the only possibility is f ■ -D which can be eliminated using the equations of 
motion [^). For spin-| baryons, there are eleven such operators 



jk 



j^JjWMi^_y,,{rj,+r,^,)^kjij^^ii^y^jj'Qi'j'k ^WM (^_y^{rij+Vj')^kjij^n'y^ii/Qi'j'k 

+6f (BBM+) str(W+) + 6f [BM+B] str(W+) + bf^^ (BB) str(W+A^+) 
+b^^ (BB) str(W+)str(A<+) + b^^ {BBW+) str(7W+) + (BW+B) str(A^+) 

(27) 

involving W+® Additionally there are eleven analogous operators involving the opera- 
tor combination W_ (S> Ai-, however these terms do not contribute to the baryon masses at 
this order. Similarly the eight operators containing v ■ >V_ do not contribute to baryon 
masses to 0{e'^). 

For the spin-| baryons, there are six operators at 0{mQ a) that contribute to the baryon 
masses. These terms are 



Av 



+^3^^ {T,M+T,) str(W+) + tY'' (r,T,) str(W+A^+) 

str(>V+)str(A^+) + £Y^' (r^W+r^) str(A^ + 



(28) 



Additionally there are six analogous terms involving the operator combination >V_ ^ Ai- 
and four terms containing v ■ A® >V_, but these operators do not contribute to the baryon 
masses at this order. 

Next we must assess the contribution from the different classes of operators in C^^\ 
Class one operators are quark bilinears that do not break chiral symmetry in the valence 
and sea sectors. These operators, however, break the full graded chiral symmetry down to 
the chiral symmetry group G of the mixed lattice action. To describe such operators in the 
effective theory, we rewrite the class one operators of the Symanzik Lagrangian in a different 
form. Instead of describing the bilinears in terms of wq and wq matrices, we use the linear 
combinations Wq + wq and Wq — wq. The former combination is just the identity and terms 
in the effective theory that stem from it are trivial to construct because they are chirally 
invariant. Those with wq — wq have the flavor symmetry of the mixed action. Let us write 
out one such term from C^^\ 

Qlp'iwQ - wq)Q = QLlp'iwQ - wq)Ql + Qnlp'iwQ - wq)Qr (29) 
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To make these terms SU{4\2) chirally invariant, we introduce two spurions Oi and O2 
transforming as 

01 LOiL^ 

02 ^ R02R\ (30) 

that will be assigned the values Oi = O2 = wq — wq. In the effective field theory, the 
operators O = C.'^OiC,, and ^^2'^^ both transform as UOU^ . Instead of using these operators, 
we work with linear combinations of definite parity, namely 



(31) 



Thus class one operators get mapped into the effective theory as terms contained in the 
Lagrangian 



bo [BB) + 6f {BBO+) + [BO+B) + 6f [BB) str {O^ 



-to {T,T,) - tf {T,0+%) - {T,T,) str {O, 



(32) 



When the action is unmixed, the operator is proportional to the identity matrix and 
hence there is only one independent operator for the ;B-field and one for the T^-field in 
Eq. dni. 

All class two operators have an insertion of the quark mass matrix ttlq. Thus these 
operators are at least of 0{mQ a^) = 0{e^) in our power counting and can be neglected to 
the order we are working. Class three operators are four-quark operators that do not break 
chiral symmetry in the valence and sea sectors. They do break the SU{4\2) chiral symmetry 
down to the chiral symmetry G of the mixed action. The requisite spurions for these class 
three operators are Oi<^Oi, 0\®02-, 02®0\, and 02®02- Hence operators in the effective 
theory will involve the products (9+ ® 0+ and O- ® 0-. The latter do not contribute to 
baryon masses to the order we work, while the former baryon operators are contained in the 
Lagrangian 



+6f ^ {BBO+) str(C+) + (BO+B) str(0+) + b^^ (BB) str(C+0+) 
+6^^ {BB) str(C»+)str(0+ 



for the spin- 1 fields, and 



(33) 



+1^"" {T^O+T,) str(0+) + ^ (r^r^) str(0+0+) + ^ (T^T^) str(0+)str(OH 



(34) 
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for the spin- 1 fields. When considering unmixed actions, all of the operators in Eqs. fl33|) 
and become redundant compared to those in Eq. because 0+ is proportional to 
the identity. 

Class four operators are four-quark operators that break chiral symmetry in the Wilson 
sector of the theory. Such terms involve the flavor structure wq ® wq, which must be 
promoted to various spurious. In the Symanzik Lagrangian, for example, we have the 
operator 

{QwqQY = {Ql^qQrY + {Ql^qQr) {QrWqQl) 

2 

+ {QrWqQl) {QlWqQr) + {QrWqQl) , (35) 
and so we require the spurious [l^ 

Bi®B2 LBiR^ ® LB2R^ 

b\®bI rbIl^^rbIl^ 

Cl ® C2 ^ RCiL^ ® LC2R^ 

cl ® cl LCIr^ ® RCIl\ (36) 

that will ultimately be given the values Bi = B2 = Ci = C2 = aAgfj^wg, and similarly 

for their Hermitian conjug ates. Now the operators O = ^^5l,2^^ ^bI^^, ^C^^, ^t^t^t^ 

^^^2^^, and ^Cl^ all transform as UOW under their respective spurion transformations. 
Moreover when assigned constant values for their spurious, the operators involving C's and 
i?2's become indistinguishable from those involving Bi and Finally instead of working 
with the operators C,wqC, and C,'^wq^\ we work with linear combinations that are parity even 
and odd, which are the W± operators, respectively, that were introduced previously. Thus 
our spurion analysis shows terms in the effective theory will involve products of Wilson 
operators W+ ® W+, and W_ ® W_.^ The latter products do not contribute to the baryon 
masses at this order. Thus for the spin-i baryons, we have the eight terms from class four 
operators 



(W+W+)''''' B'^''' + bYs''^' (W+W+)"' B''^'' 



^QCD 

+6f {BBW+) str(W+) + {BW+B) str(W+) + b^ {BB) str(W+W4 
+bf (BB) str(W+)str(W+) 



(37) 



6 



As spurions, these have the same effect as squares of the 0(a) spurions. Thus we need not consider such 
higher-order effects from lower order terms. 



11 



while in the spin-| sector, we have five terms 



A 



QCD 



)'Tk3i\M 



+tY [r^W+r^] str(>V+) + tY {T^T^) str(>V+>V+) + {r^T^) str(>V+)str(W+) 



(38) 



Class five operators break the 0(4) rotational symmetry of Euclidean space. The lowest- 
order hypercubic invariants that are parity even are: v^v^v^v^^ v^v^S^S^, and S^S^S^S^. 
The latter two are both redundant because, suspending the Einstein summation convention, 
we have S^S^ = \{Sfj^^ + v^v^). There is an additional hypercubic invariant for the spin-| 
fields since they carry a vector index. Furthermore, these class five operators have the chiral 
symmetry of the group G of the mixed action, not the full graded chiral symmetry of C^^\ 
As with class one operators, class five operators in the effective theory require the insertion 
of the operator 0+. Thus the 0(4) breaking operators are 



C 



B) + h\ [B 



BO+) + [B 



O+B) 



+63 (Bv^v^Vf^v^B) str (C+) - {T^v^v^v^v^%) - tl (T^v^v^v^v^O+T^) 
-tl (T^v^v^v^v^T^) str (0+) - [T^v^v^T^,] - tl [T ^v^v^^O+T^) 



-tl [r^v^v^T^) str (0+) 



(39) 



The mixing of classes of operators is only possible through linear superpositions. This is 
accounted for in the effective theory by linear superpositions of effective operators and thus 
need not be addressed separately. This completes the (nearly) exhaustive listing of 0{e'^) 
PQxPT terms for the mixed lattice action. Despite the large number of operators (over one 
hundred), many terms do not contribute to baryon masses at tree level and have already 
been omitted. Additionally many of the contributing terms are not linearly independent at 
leading order. In the following section we calculate the masses of the nucleons and deltas to 
(9(a^), and thereby determine the number of new free parameters that are required. 



III. BARYON MASSES TO 0{a^) 

In this section we calculate the masses of the nucleons and deltas in SU{A\2). Baryon 
masses in SU{2) are presented in Appendix 1X1 



A. Nucleon mass 

In Euclidean lattice field theory, the baryon self energy E(7^p^) is no longer rotationally 
invariant; it is a hypercubic invariant function of the baryon momentum and gamma matri- 
ces. With this functional dependence in mind, we can now calculate the mass of the nucleon. 
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FIG. 1: Loop diagrams contributing to the a-dependence of the nucleon mass at 0(6^). Mesons 
are denoted by a dashed Une, flavor neutrals (hairpins) by a crossed dashed line, and a thin solid 
line denotes an octet baryon. The square denotes an 0{a) vertex. 

The nucleon mass in the combined lattice spacing and chiral expansion can be written as 

= Mo ifi) - M^') (/i) - M^p ifi) - Mj;) (/i) + . . . (40) 

Here, Mq (/x) is the renormalized nucleon mass in the continuum and chiral limits, which 
is independent of tuq. is the contribution to the mass of order e", and /i is the 

renormalization scale. 

At order e"^ in the combined lattice spacing and chiral expansion, we have contributions 
at tree level from the 0{mQ) and 0{a) operators in Eq. (|22j) . These contributions to the 
nucleon mass read 

M^^^ = 2{aM + f3M)mu + rrij + 2aAQ^£, [{aw + + 2aw Wg] • (41) 

From this expression, we see that the 0{a) mass corrections vanish only if both the valence 
and sea quarks are GW. This is in contrast to the meson sector where only the valence 
quarks need be GW to ensure the vanishing of linear a contributions. At 0{e^) there are 
contributions from loop graphs to M^"*. These have been determined for SU{4:\2) in [i^l 
and we do not duplicate the expressions here because the only modification necessary is to 
include the a-dependence of the loop meson masses via Eq. (fTTf).'' 

At 0{e'^), we have the usual continuum contributions from tree-level 0{mQ) operators 
and from loop diagrams. These have been detailed for S'f/(4|2) in js^ and we do not 
duplicate these lengthy expressions here. The only modification necessary at finite lattice 
spacing is the inclusion of the a-dependence of the loop meson masses from Eq. (|17p. The 



^ Additionally there are local operators that contribute at this order. The form of these operators, however, 
is the same as those in Eq. (|22|) but multiplied by a factor of A/A;^. We can thus trivially include the 
effect of these operators by promoting the LECs to arbitrary polynomial functions of A/A^ expanded out 
to the appropriate order. We do not spell this out explicitly since the determination of these polynomial 
coefficients requires the ability to vary A. 
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remaining finite lattice spacing corrections arise at tree level from the 0{mQ a) operators in 
Eq. and the 0{a?) operators in Eqs. (jHH), (|H7j). and (jSHI)- Additionally there are loop 
diagrams arising from the operators in Eq. ()22|) . where either the Wilson operator in Eq. ()24|) 
is expanded to second order or is inserted on an internal baryon line. These diagrams are 
shown in Fig. ^ Lastly there are wavefunction renormalization corrections that are linear 
in a. The contribution to the nucleon mass from these operators and loops appears in M^^^ 
as a correction of the form 



5M 



(4) 



N 



A 



X 



-4 — — — < {aw + M [{wv + Ws)C{mju, /i) + WyC{m^^ , , /i)] 



A2 

X 



+ 2aw Ws [2C{mjj, /i) + /:(m^^ , m^^^, , n)] + Af^ 



+ gUXB'^Jim^. A, /i) + B'yj{m.y,, A, /i)] 



N 



Wk\cD E + E'w^ + E"ws + E'" 

+ u{p)v^v^v^v^ u{p) {F + F'w^ + F"ws) 

The non-analytic functions appearing in this expression are defined by 



(42) 



C{m^,n) = mAog — 



5 - . 5"^ - 7111 + ie 



5 + . 5"^ - m'i + ie 



J{m^, m^/ , (5, /i) = 7^00/ [J{m^, S, /i) , J{m^' 



(43) 

(44) 

(45) 
(46) 



and one should keep in mind the a-dependence of the meson masses. The parameters A-F 
are replacements for particular combinations of LECs.^ The meson loop coefficients Af^, 



® For clarity in comparing with the operators written down in Sec. Ill CI these combinations are: A 
2bY^ + 2bY^' + bf^' + 6f + 2bf^', B = 2b^'' + 26f *^ C = 2b%^' + 2bfi^' , D = 2b 

2b^'^ + 2b^'^ + 2bf^ + 46f E'^b'^ + b^ -r U3 



WM ^ 2bY^, 



E = bo-b? -bf -2bf + b?'^ + bf'^ + b^<^ 
bf + 2bf + 2b^ ~ Abf^ ~ Ab<i°, E" - 2bY + 468^ + 46^ - Abf^ ~ 46^^, E'" = 2bf + 2bf + ^b'^^ + Sb^'^, 
F = bl-bl-bl- 2bl, F' = 2bl + 2bl, and F" = 46^. 
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, and B^^ are given by 



-2g\aw + QaQiPw 



N 



B 



N 
ju 



Wy{aw + Pw + 1w) + 'iwsiaw - crw), 

Wv (aw + + \ + Ws (^aw - 2aw + ^Iw 



(47) 



Despite the complicated form of Eq. ()42|) . there are at most only four free parameters 
depending on the mixed action considered. The reason for the complicated form is that 
the above equation shows the interplay between all possible theories with mixed actions. 
Precisely which operators are present or, on the other hand, can be eliminated is perhaps 
of academic interest because at this order many of the individual contributions cannot be 
resolved from lattice data. For example, in a theory with Wilson valence and sea quarks, 
the contribution A + B cannot be further resolved. In Table ITTl we summarize the number of 
free parameters introduced by 0{a), 0{mQ a), and O(a^) operators in the baryon sector for 
various partially quenched theories. Of course, the simplest situation occurs in a theory with 
GW valence and sea quarks. Here we see that there are O(a^) corrections to the nucleon 
self energy of two forms, one of which violates 0(4) symmetry. 

Let us investigate the violation of 0(4) rotational invariance in detail. For an on-shell 
nucleon at rest, the 0(4) breaking operator merely contributes an additive shift to the 
mass. This correction with F coefficients becomes indistinguishable from the corrections 
which respect 0(4) [these have E coefficients in Eq. ()42|) ]. For any particle with non- 
relativistic three-momentum p, however, we can see the effects of 0(4) violation in the 
dispersion relation. To write the dispersion relation compactly, let us assume both valence 
and sea quarks are GW. The same form of the dispersion relation holds for actions involving 
Wilson quarks with the addition of various 0{a) contributions. First let us ignore the 
contributions from 0(4) breaking operators. The dispersion relation has the familiar form 



M{a) + 



2M(a) 



(48) 



where M(a) = M(0) — a^Aqfj^E, and M(0) denotes the particle mass in the continuum 
limit (which remains quark mass dependent). 

To add the 0(4) breaking terms to the dispersion relation, we can use a non-relativistic 
expansion for the particle four-velocity or utilize reparameterization invariance to deduce the 
corrections to the dispersion relation. We relegate this discussion to Appendix where we 
address particles with spin less than two. Using the result in the Appendix for the spin-1/2 
nucleon, the relation between the particle energy and momentum is then given by 



= MiO)-a'A-'QCDiE + F) 



2M(0) 



1 + 



M(0) 



{E + 5F) 



(49) 



Here we have written the non-relativistic expansion to second order and have retained the 
leading lattice-spacing correction for each term. Notice for a particle at rest, we recover 
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the simple result of an shift. The kinetic energy receives corrections at (9(a^), these are 
due to the shift in mass as well as the 0(4) contribution in Eq. ()B2|) . If it were not 
for the latter contribution, the dispersion relation would retain the form in Eq. ()48|) with 
M(a) = M {Q) — a? Kq(j j~,{E + F) . Additionally the relativistic correction to the kinetic energy 
will be modified with a coefficient proportional to and a different linear combination of 
E and F . Also the energy is only a cubic invariant function of the spatial momenta pi. 
Because four factors of the momentum are required, the leading 0(3) breaking term in the 
energy occurs at (9(a^M~'^) with a coefficient — F. 

B. Delta mass 

For the masses of spin-| baryons, there is an additional feature on a hypercube due to 
the fact that the Rarita-Schwinger fields themselves carry a vector index. The one-particle 
irreducible self energy for a spin-| field on a hypercube has the form 

u^{p) S^^ Uy{p) = u^{p) (S + S^^) M^(p), (50) 

where S and E^^ are both hypercubic invariant functions. With this in mind, we can now 
calculate C(a), 0{mQa) and 0{a?) contributions to the delta mass. The mass of the delta 
in the combined lattice spacing and chiral expansion can be written as 

Ma = Mo (/x) + A(/x) + ^ (/i) + Mf> (/x) + uf (/i) + . . . (51) 

Here, Mq (/i) is the renormalized baryon mass in the continuum and chiral limits and A(/i) is 
the renormalized mass splitting between the spin-| and spin-| baryons in the continuum and 

chiral limits. Both of these quantities are independent of a and ttlq. is the contribution 
of order e", and /i is the renormalization scale. 

At order e"^ in the expansion, we have contributions to the delta mass at tree level from 
the Oimq) and 0{a) operators in Eq. (|^. These lead to 

M'f' = 27m rriu - rrij + 2aAQ^,j^ [jw Wv - '2a w Ws] . (52) 

Again in contrast with the mesons, we note that both valence and sea quarks must be GW 
for the 0{a) delta mass corrections to vanish. At 0{e^) there are contributions from loop 
graphs to M^-*. These have been determined for 5'?7(4|2) in jH^] and we do not duplicate 
the expressions here, because the only modification necessary is to include the a-dependence 
of the loop meson masses via Eq. (fT7j) . 

Finally at 0{e'^) there are the usual contributions from tree- level O(mQ) operators and 
from loop diagrams. For the delta, these have been detailed for SU{A\2) in j^. The 
remaining finite lattice-spacing corrections arise at tree level from the 0{mQ a) operators in 
Eq. (|2HI) and the 0{a?) operators in Eqs. (1221), (EHI), and dSHI)- Additionally there are loop 
contributions stemming from the operators in Eq. (j22|) . These loop diagrams are depicted 
in Fig. 121 Lastly we must include wavefunction renormalization corrections that are linear 
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FIG. 2: Loop diagrams contributing to the a-dependence of the delta mass. Mesons are denoted 
by a dashed line, flavor neutrals (hairpins) by a crossed dashed line, and a thick solid line denotes 
an decuplet baryon. The square denotes an 0{a) vertex. 



in a. The contribution to the delta mass from these operators and loops yields 

AWyTTlu + EWsTTlu + CWyTTlj + DlUsTrij 



5^(4) ^ ^KCD 



A 



X 

aA2 



QCD 



A2 
X 



— 2(7 w Ws ['^J^i'lTT'jj, /^) + >C(m 



^ 26 



- \91n [B^Ji^n. - A, /.) + Bf^J{m,u, -A, /.)] | 



WaIcd E + E'wy + E"ws + S'" 

+ u{p)uVnV^v^v^ u,,{p) (f + F'ly^ + F" 
+ Uf,{p)v^Vf, u^{p) (g + G"w^ + C'ty^^ 

The meson loop coefficients and B^ are given by 

B^ = ^^,(011/ + /^w + 7w) + 2wsiaw - o'w), 



JU 



— -['Wy{5aw + + 67^) + Ws{aw + 4/3iy + I2aw - I2aw)]- 



(53) 



The complicated nature of this expression results from considering general theories with 
mixed actions. The formula retains book-keeping factors that indicate which operators 
contribute in particular mixed action theories. Notice, for example, that for a theory with 
Wilson valence and sea quarks, the combination A+B cannot be further resolved from lattice 
data. Here the parameters A-G depend on particular hnear combinations of the LECs in the 
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TABLE II: Summary of the new free parameters entering particle masses in chiral perturbation 
theory for non-zero lattice spacing. Listed for PQ^PT and xPT are the number of linearly inde- 
pendent operators contributing to the masses of the pion {ttIt^), nucleon {Mj\f), and delta (Ma). 
These operators are classified as 0{a), 0{mQ a), or O(a^). For each theory the quarks are grouped 
by species and for partially quenched theories are grouped by pairs of valence and sea quark species. 
The species are labeled W for Wilson, and GW for Ginsparg- Wilson. 



PQxPT 










Mn 






Ma 




(v) is) 


0{a) 


0{mQ a) 




0{a) 


0{mQ a) 




0{a) 


0{mQ a) 




WW 
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3 


W GW 
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1 
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SU{2) xPT 


0{a) 


0{mQ a) 




0{a) 


0{mQ a) 




0{a) 


0{mQ a) 




W 


1 


1 


1 


1 


1 


2 


1 
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3 


GW 

















2 
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spin-| sector.^ The number of free parameters introduced for Wilson valence and sea quarks 
from 0{a), 0{mQ a), and 0{a'^) operators is six. Table ITTl lists the number of free parameters 
introduced from finite lattice spacing effects for delta mass in mixed action theories. The 
simplest expression results from an unmixed action of GW quarks. Corrections of course 
start at O(a^), but are of three distinct forms given the possible hypercubic invariants. 

As we know from determining the nucleon mass above, the 0(4) violating terms propor- 
tional to v^v^v^v^ in Eq. with coefficients F give rise to corrections to the dispersion 
relation. There are additional corrections with coefficients G that are tied to the Rarita- 
Schwinger spinors in the form u^{j))v^v^ u^{p). Because of the constraint u^{p) = 0, these 
new terms do not contribute for a particle at rest. To consider a particle moving with non- 
relativistic three momentum p, we must use explicit forms for the Rarita-Schwinger spinors 
to obtain the modification to the energy-momentum relation, see Appendix |B] For GW 
valence and sea quarks, we find a similar dispersion relation given in Eq. ()49|) but with an 
additional term 

4 = M(0) + aV,UE + F) + ^ [l - f^(E - 3F) 



■^[|A|pi-(2|A|-3)Pa. 



where A = ±3/2, ±1/2 refer to the delta spin states quantized with respect to an axis, p\ 
is the momentum along that axis, and p± is the momentum transverse to that axis. We 
have retained only the leading term in the non-relativistic expansion and its correction. 



9 For clarity in comparing with Sec HTCl these combinations of LECs are: A = 21^"^ + tz^*^, B = 2t^^^, 
C = 2tf, D = 2tY^' + 4tf ^^ £; = io - if - 2tf + if « + if « + 2t^^ + 21^^ + Atf<^, E' = 
tY +tY + 2if - 4if^», E" = 2tf + Atf + 4tf - U^o^ e'" = 2tf + 8tf F ^ - t\ - 2q, F' - 21^1, 
F" ^Atl. G^tl-tJ- 2tl, G' = 211, and G" = 4iJ, 
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Unlike the nucleoli, the delta dispersion relation is more sensitive to 0(3) violation which 
enters at 0{a^M~^) as opposed to 0{a^M~*). The splitting of spin states is proportional 
to |A| because one has invariance under rotations of the quantization axis by vr. For an 
unpolarized state, however, the 0(3) violation in Eq. averages out. 

IV. SUMMARY 

Above we have extended heavy baryon chiral perturbation theory for the Wilson action to 
O(a^). In our power counting, we have included all operators that are at least of 0{e*), this 
includes operators of Oimqa) and O(a^). We have considered partially quenched SU{A\2) 
with a mixed action, as well as SU{2) in Appendix^ To this order in the combined lattice 
spacing and chiral expansion, we saw the necessity for the introduction of a large number 
of new operators. For the case of a partially quenched theory with a mixed lattice action, 
there are well over one hundred operators in the free Lagrangian containing the spin-i and 
spin- 1 baryons. 

Despite the introduction of a large number of terms with undetermined LECs, observables 
calculated to O(a^) depend on still relatively few independent parameters. This is because 
many of the operators introduced for mixed lattice actions act differently depending on 
the number of sea quarks contained in particular states of the multiplet. Each level of the 
multiplet must be treated differently in the effective theory because there are no symmetry 
transformations between the valence and sea sectors of the theory. This distinction is impor- 
tant if the level-dependent operators act on baryons within loops and such a situation only 
occurs at higher orders in the chiral expansion. Moreover, most of the remaining operators 
for baryons which consist of only valence quarks are not independent at leading order and 
hence cannot have their LECs disentangled from lattice data. To make these points clear, 
we have calculated the finite lattice-spacing corrections to baryon masses in various theories. 
In Table IHl we summarize the number of free parameters due to lattice spacing artifacts 
that enter into expressions for baryon masses and can be determined from lattice data. The 
table is separated for the pion, nucleon, and delta into contributions from independent 0{a), 
0{mQ a) and 0{a?) operators. 

There are some further points to observe about the baryon sector at O{o?'). While 
there are considerably more operators in the partially quenched theories for mixed actions 
as compared to the unquenched theories, the number of independent parameters entering 
into the determination of baryon masses from lattice data is the same for 0{a) and O{o?) 
operators; there is only one additional 0{mQ a) parameter which stems from the different 
sea quark mass. Nonetheless, the computational benefits of partially quenching are not 
hindered by the explosion in the number of new baryon operators because most are not 
independent to 0{e^). 

On the other hand, based on studies in the meson sector [l^, one would expect the 
number of free parameters to be reduced when considering theories with Ginsparg- Wilson 
valence quarks and Wilson sea quarks compared to the Wilson action, cf. Table |TT| While 
there is a clear and sizable reduction in the number of baryon operators for actions employing 
Ginsparg- Wilson valence quarks and Wilson sea quarks compared to the unmixed Wilson 
action, there is no reduction in the number of free parameters involved in the baryon masses. 
The only reduction in the number of free parameters comes about in simulations employing 
Ginsparg- Wilson quarks in both the valence and sea sectors. Another notable difference for 
heavy baryon fields is that operators which reduce the 0(4) symmetry to the hypercubic 
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group are present in the effective theory. Such operators were suppressed in the pseudoscalar 
meson sector because the absence of a large mass scale implies the pseudoscalar momenta 
are small. For heavy baryon fields, however, there is a large mass scale and the momenta 
are automatically large, thus such operators are not suppressed when multiplied by powers 
of the lattice spacing. Similar operators must then enter at 0{e'^) for other heavy particles. 
We have spelled out the form of 0(4) breaking corrections for particles of spin less than two 
in Appendix FBI 

In this work we have seen that near the continuum limit, the spin-| baryon self energy 
has the behavior^" 

M(a) = M(0) + a a + P anig log + 7 a + 5 + e u{p) Vf_,v^v^v^ u{p) + . . . , (55) 

where the omitted terms are of order and higher. In writing the above expression, we have 
assumed for simplicity that one fine tunes the quark mass so that the pion mass vanishes 
in the chiral limit. In this case there are only polynomial corrections in a. The influence 
of the 0(4) breaking term on the baryon energy was determined in Sec. IIII Al For spin-| 
baryons, the behavior of the self energy has the same form, however there is an additional 
contribution to the self energy of the form 

Uf,ip) S^^ Uf,{p) =uja^ u^{p) v^v^ u^{p) + (56) 

and we have detailed the form of this correction in Sec. IIII Bl A similar such additional term 
is present for vector mesons as well. 

Knowledge of the quark mass dependence of baryon observables is crucial to perform the 
chiral extrapolation of lattice data, extract physical LECs, and make predictions for QCD. 
Artifacts of approximating spacetime on a discrete lattice make the chiral extrapolation 
more challenging because of the introduction of additional error. With xPT and PQxPT 
formulated for the Symanzik action, one can parametrize the dependence on the lattice 
spacing and considerably reduce the uncertainty surrounding lattice artifacts. The formalism 
set up here is straightforward to extend to other Wilson-type fermion actions, twisted mass 
QCD [5^, for example. On the other hand, due to the proliferation of operators, it is 
doubtful that partially quenched baryon staggered xPT would provide much analytic insight 
at 0{a^). 
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APPENDIX A: BARYON MASSES TO 0{a'^) IN SU{2) xPT 

In this Appendix, we detail the case of finite lattice-spacing corrections to baryon masses 
in SU{2) xPT. First we write down all finite lattice spacing terms that arise to 0{e'^) from 



For GW quarks, the coefficients a, (3, and 7 are identically zero. 
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the Symanzik Lagrangian. Next we determine the corrections from these operators to the 
masses of the nucleon and delta. 

The form of the Symanzik Lagrangian for SU{2) flavor is mainly the same as in the 
main text. Of course there is no separation of the theory into valence and sea sectors. The 
mass matrix for two light flavors in the isospin limit is given by mq = diag(m„,m„), while 
the Wilson matrix appears as wq = diag(w^,w„). The subscript v now has no particular 
significance in SU{2) xPT other than to maintain the definition of Wy used in above. 

At zero lattice spacing and zero quark mass, two flavor QCD has an SU{2)l ® SU {2)r ® 
U{l)v chiral symmetry that is broken down to SU{2)y ® U{l)v- xPT is the low-energy 
effective theory that emerges from perturbing about the physical vacuum. The pseudo- 
Goldstone bosons can be described by a Lagrangian that takes into account the two sources 
of chiral symmetry breaking and is given in Eq. ()14|) . with the exception that $ is an SU{2) 
matrix containing just the familiar pions, and the str's are now tr's. Their masses are given 
in Eq. (ITTIl . 

In SU (2) xPT, the delta baryons are contained in the flavor tensor TJf'^ which is embedded 
in the PQxPT tensor T^^'^ simply as T^^^ = TJf'', when all indices are restricted to 1-2. 
Consequently the form of the xPT Lagrangian for the delta fields has a form very similar to 
that of PQxPT. The nucleons, however, are conveniently described by an SU (2) doublet 

These states are contained in the PQxPT flavor tensor B^^'' as |4 

V6 

when all of the indices are restricted to 1-2. Consequently operators involving the nucleons 
will appear differently in SU{2). 

To 0{e'^) the free Lagrangian for the nucleons and deltas reads 

£ = iNvDN - 2aNmi{M+) - 2aJ^mi{W+) 

+iT^v ■ DT^ + AT^T^ - 2aT,,TMM+) - 2a^T^r^tr(W+) . (A3) 

The LECs for the nucleons and deltas in S'[/(4|2) are related by the matching equations 

= -]pM + ^M, (A4) 

and completely analogous equations that relate cx^, and to the PQxPT parameters a^y, 
I^Wi <^w-, IWi and 'ow- The interaction Lagrangian of SU{2) appears as 

C = 2gANS ■ AN + 2g^N (T^A^N + NA^T^) - 2g^J'^S ■ AT^. (A5) 

The relation of the parameters appearing in the PQxPT Lagrangian in the main text to 
those in xPT can be found from matching. One finds ji^]: Qa = ^c^m — ^Pm, San = —C, 
and gAA = 'H. 
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The first lattice spacing corrections at 0{e^) arise from tlie 0{mQ a) operators. In de- 
generate SU{2), we liave only two such baryon operators 



Av 



^wM NNtiiW+)ti{M+) -t^^^T^T^tT{W+)tiiM+ 



(A6) 



The remaining terms are of order O(a^) and appear in the Lagrangian 
1 



A 



QCD 



NN ii{W+)ti{W+) - t^T/,T^tr(W+)tr(W+) 
n NN + n"Nv^v^v^v^N - t T^T^ - f T^v^v^v^v^T^ - f T^v^v^T^ [A7) 



The coefficients of operators in the SU{2) theory are contained in SU{A\2), but because 
SU{4:\2) contains considerably more operators, the matching yields algebraically cumber- 
some expressions that relate the LECs. 

Having written down all of the relevant operators for the nucleon and delta masses, we 
now calculate their lattice spacing corrections. We do not write down the loop contributions 
from a-independent operators because the only modification necessary is to include the 
lattice spacing dependence of the meson masses that appears in Eq. (jl7|) . For the mass of 
the nucleon, at O(e^) we find 



M; 



(2) 



N 



4cr ruu + AaAgc^a^, Wy, 



(A8) 



The finite lattice spacing corrections to M^^"* read 



SM 



(4) 



N 



4w. 



aA 



QCD 



X 



n^'^rriu - 3a^C{m^, /i) + ^^^(cr^ - a^) [j{m^, A, /i) + ml] 



+a Aqcd \n + 4:71 Wy + rfu{p)v^v^v^v^ u{p)j ■ 
For the deltas, similar calculations yield the result 

^ = -4CTm„ - AaAQcoO'wWy. 
At 0{e^) the remaining finite lattice spacing corrections to delta mass are 



(A9) 



(AlO) 



6M 



(4) 



A2 

X 



A^ t^^TJiu + ?>ayjC{m^, /i) + 5'AAf - (yw)J{m^, -A, /i) 



+a^A^c'D [t + 4t Wy + t''u„{p)v^v^Vf,Vf, u^{p) + t"u^,{p)v^v^ u^{p) ) . (All) 



APPENDIX B: 0(4) BREAKING OPERATORS 

Above we addressed 0{a^) corrections in the baryon sector. At this order we saw that 
0(4) breaking operators are present in the chiral effective theory. In this Appendix, we 
detail the form of these corrections in heavy particle effective theories for particles of spin 
less than two. 
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We consider the simplest case first, that of a heavy scalar field 0. At (9(a^), we have an 
0(4) breaking operator of the form 

£ = a^co (j)^ v^v^v^v^ 0. (Bl) 

We can either treat the four-velocity as the physical velocity of the (f) field, or consider the 
four- velocity to be fixed in the rest frame, for example. Clearly when the particle is at rest, 
both descriptions are identical; and with = (0,0,0,2), we obtain a trivial dependent 
shift to the mass. Away from rest the effects of 0(4) breaking are no longer trivial. 

First let us consider to be the physical velocity of the 0. In terms of the particle's 
spatial momentum p, we have = {'yp/M, 27), where M is the mass of the and 7 is the 
Lorentz-Fitzgerald contraction factor. In a non-relativistic expansion, we can express the 
particle's four-velocity in terms of the spatial momentum. This enables us to write 

J2v,V,V,V, = l + 2^ + —iJ2p^P^P^P^+P']+0{M-'). (B2) 
ti=l \i=l ) 

Thus the operator in Eq. (jBljl . besides shifting the mass by a^co, modifies the particle's 
kinetic energy as well as the relativistic correction to the kinetic energy. The energy of the 
is now a cubic invariant function of the spatial momentum. The breaking of 0(3) rotational 
symmetry for a particle's energy-momentum relation is, however, small for non-relativistic 
momenta. 

Extending the above analysis to non-zero spin involves boosting the various spin wave- 
functions. On the other hand, if one treats the four-velocity as fixed and one can derive 
a string of fixed coefficient operators using reparameterization invariance |5^. For a fixed 
four- velocity, we write the heavy scalar momentum P as = Mf^+p^. This decomposition 
is arbitrary as one can make a reparameterization of the form 

Accordingly the field 0(x) changes to e*'^''^^0(x), and to maintain the normalization of the 
new four- velocity, we require v ■ q = — g^/(2M). The physics is unchanged by this transfor- 
mation and this is referred to as reparameterization invariance. In order for the Lagrangian 
in Eq. ()B1|) to be reparameterization invariant (RPI), there must be additional operators 
with fixed coefficients. These can all be written compactly in the manifestly RPI form 

--^^o.t(„„.^)(...^)(„„.^)(.„.^).. ,B3) 

Expanding out the terms of Eq. ()B3|) and acting between states of residual momentum 
taken in the non-relativistic limit, we recover the result of Eq. flB2|) . 

For a heavy spin- 1/2 particle tp, we have the 0(4) breaking Lagrangian 

C = a^ci/2ip v^v^Vf^Vf, ip, (B4) 

which is not RPI. Unlike the scalar case, we have a further constraint on the field due 
to the spin degrees of freedom. To maintain the condition —i'^ijj = ijj, the field ?/'(a;) must 



23 



change under the reparameterization in Eq. ()B3|) . To order M ^ , we have reparameterization 
invariant spinor \l/ given by js^ 

(B5) 

where -D;x,± = + v^v^D^,. Since the bihnear = + 0{M~^), the RPI form of 
Eq. dBD) is 



M J M J y"^ M J m ^ 

and hence the correction to the energy to 0{M~'^) has the same form as in Eq. ()B2|) . If 
we were to write the covariant spinor to fourth order, we could ascertain the coefficient of 
the p'^/M^ term in Eq. (jB2l) for a spin-1/2 particle. The 0(3) breaking term, however, is 
unchanged. 

The correction for a heavy vector particle, the rho meson for example, is deduced 
similarly. The O(a^) Lagrangian contains two terms 

C = a^ci pI v^v^v^v^ py + a^ci p\ v^v^ p^, (B7) 

and hence there is a new type of contribution to the energy proportional to e* (p) v^^v^ ^fj,{p), 
where e^(p) is the rho's polarization vector. This correction vanishes for a rho at rest. The 
terms in Eq. ()B7|) are not RPI. To deduce the RPI form of the Lagrangian, we must now 
also maintain the constraint v^p^ = 0. We can build the RPI Lagrangian from RPI vector 
field Rf,. To C(M-2), we have 

= (s,. - + ^) P., (B8) 

and hence the RPI form of the Lagrangian in Eq. ()B7|) can be deduced from combinations of 
RjiTRfj,. Because we are addressing the corrections at tree level, the states are on-shell and 
the additional terms in Eq. ()B8|) vanish. In loops the off-shell degrees of freedom propagate 
and the additional terms are necessary. The on-shell RPI Lagrangian is thus given by 

pI + ^) (^,, + ^) p„. (B9) 

We can use the above Lagrangian to determine the O(a^) corrections to the vector parti- 
cle's dispersion relation. Let us denote M(a) as the particle mass including lattice discretiza- 
tion contributions. One such contribution is a^ci that arises from the above Lagrangian. 
There are many others, however, that we have not spelled out. Modification of the disper- 
sion relation only arises from the terms in Eq. ()B9|) . Utilizing the non-relativistic expansion 
for the momentum p and explicit forms of the polarization vectors, we arrive at the energy 
of a spin-one particle 
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where M(0) is the rho mass in the continuum hmit, A = ±1, is the spin along an axis, p\ 
is the momentum along that axis, and p± is the momentum transverse to that axis. For a 
vector particle of given polarization, the energy- momentum relation breaks 0(3) rotational 
invariance at 0(a^M~^). Notice for an unpolarized state, the rotational symmetry breaking 
from this term is averaged out and then does not enter until (9(a^M"^) as with scalar and 
spin-1/2 particles. 

Finally we address the corrections for the case of a spin-3/2 particle represented by a 
Rarita-Schwinger field ipfj,. At O(a^) terms that break 0(4) appear in the Lagrangian as 

C = a^C3/2 v^v^v^Vf, ipu + a%/2 ip^ v^v^ ip^. (Bll) 



To deduce the corrections to the energy, we cast this Lagrangian into its RPI form. Under a 
reparameterization, one must maintain the conditions —i'^ipfi = V'm '^fj.'^n = 0- The RPI 
field is a combination of terms from RPI invariant spinor and vector fields and is given 
by 



(B12) 



up to 0{M^^). Half of the terms in the RPI spin-3/2 field vanish on shell. The remaining 
terms are identical to those for a spin-1/2 particle. Thus for an on-shell particle, we have 
^At^/x ~ '^fi'^f^ + 0{M~^). Thus the on-shell RPI Lagrangian corresponding to Eq. ()B11|) is 
given by 

+ ^) + ^) -A.- (B13) 

As with the spin-one case, we can evaluate the corrections to the spin-3/2 energy in the non- 
relativistic limit using the Lagrangian Eq. ()B13|) and explicit forms of the Rarita-Schwinger 
vectors. The result for the A mass is presented in the main text, see Eq. (jSH). 
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